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Abstract 

A general, system-independent formulation of the parabolic Schrodinger-Poisson equation is 
presented for a charged hard wall in the limit of complete screening by the ground state. It 
is solved numerically using iteration and asymptotic-boundary conditions. The solution gives a 
simple relation between the band bending and charge density at an interface. I further develop 
approximative analytical forms for the potential and wave function, based on properties of the 
exact solution. Specific tests of the validity of the assumptions leading to the general solution 
are made. The assumption of complete screening by the ground state is found be a limitation; 
however, the general solution still provides a fair approximate account of the potential when the 
bulk is doped. The general solution is further used in a simple model for the potential profile of an 
AIN/GaN barrier, and gives an approximation which compares well with the solution of the full 
Schrodinger-Poisson equation. 

PACS numbers: 73.20.-r,71.20.Nr,74.78.Fk,03.65.-w 
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Quasi two-dimensional electron gases (2DEGs) form at many planar interfaces and sur- 
faces where electron accumulate in inversion layers.-^ They play a central role for the op- 
eration of many devices, for instance for metal-oxide semiconductor (MOS) devices, and 
high-electron mobility transistors (HEMT). Naturally, their properties such as quantized 
levels and conduction band bending have been much studied both experimentally and 
theoretically.- - — In particular the angle-resolved photoemission spectroscopy (ARPES) char- 
acterisation of the InN surfaces has spurred recent activity.— - — Hetero junctions of highly 
polar materials, such as the III-V nitrides, 15 induces these 2DEGs at the positively charged 
interfaces. A good account of the band bending at interfaces in these materials is essential for 
band-gap engineered intersubband devices such as resonant-tunneling diodes and quantum- 
cascade lasers. Simple quantum-mechanical systems, such as the particle in box, harmonic 
oscillator, and linear potential well are instructive model useful to generate rough accounts 
of various physical phenomena described by the Schrodinger equation. In the same vein, 
the charged hard wall represent a model case for the Schrodinger-Poission (SP) equation 
describing the quantization and band bending at interfaces. 

The conduction-band edge, or potential, V and quantized levels E n at interfaces are 
usually obtained with the SP equation with mass m, dielectric constant e, 

^{z) = E^{z), -e^V(z) = p(z) , (1) 

where p(z) is the total charge-density comprised of donor, interface, and electron charge. The 
related textbook linear-potential well problem is inappropriate because it lacks an account of 
the electron screening inherit to the problem. The SP equation is usually solved iteratively; 
V is updated until it reaches self-consistency. This approach is straightforward to implement, 
but as a first line of attack to device modelling and for understanding physical trends, simple 
analytical results are also of great value. 

In this paper, a general, system-independent, formulation of the parabolic Schrodinger- 
Poisson equation for a charged hard wall is presented in the limit of complete screening of 
the interface charge by the ground-state, that is, in the quantum electrical limit .-^ It is 
solved numerically using iteration. These steps follow the earlier work of Pals,- who also 
provided an analytical approximation using the variational principle. In contrast to Pals, 
I present analytical expressions that are based on constraints from physical principles and 
exact properties obtained from the numerical solution. Furthermore, I make specific tests 
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of the general solution outside its expected range of validity. Finally, I demonstrate the 
usefulness of the analytical expressions for making simple models of the band bending in 
AIN/GaN heterostrucutes. 

A key assumption made to arrive at the model system is the infinite potential barrier or 
hard wall at z = 0. It is appropriate for flat surfaces, as the potential variation is abrupt 
and the work function is much larger than other characteristic energies; for interfaces, the 
band offset must be large. Another, is the neglect of non-parabolicity, which is an important 
effect in some semiconductors, but more-so for excited states of narrow quantum wells than 
for the ground-state of the shallow quantum wells that form at interfaces. 

The assumption of complete screening of the interface charge a by the ground state ipQ 
leads to 

p{z) = a[8(z) - \M*)\ 3 ) , (2) 

with a = m(eF — Eq)/ti for an isotropic 2DEG in zero magnetic field. This assumption is 
a serious limitation, as it is both a zero-temperature (T = 0) condition and a restriction on 
the amount of charge at the interface. For T = 0, it is valid when the Fermi level is below 
the first excited state. 

The dimensionless equation that leads to the general solution are obtained with the 
change of variables: E = jK, z = Ax, and <fi = ipy/X. Here A = (h 2 e/2ma) 1 ^ 3 and 7 = Aer/e 
defines length and an energy scales. The two first change of variables are identical to the 
textbook procedure for a linear potential well.— We get 
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and note that we need only solve this equation once. The system-specific wave function if}, 
potential V, and ground state eigenvalue E can be restored for specific values of a, m, and 
e. 

To guide the computational procedure, I first consider certain limits. For large x, \K\ » 



\U\, and the asymptote of the wave function follows <fi — > exp [— yj — Kx\ , which further 
leads to the asymptote of the potential U — > Au exp[— gx], where Au = A 2 ,/AK and g = 



-2\J—K is the decay factor. For small x, the wave function is linear to first order <p(x) 
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FIG. 1: Numerically determined wave function and potential for the system-independent formula- 
tion of the Schrodinger-Poisson equation for a charged hard wall. The full dark line in the upper 
(lower) panel gives the wave function (potential). The dashed dark [blue] lines show the asymptotic 
curves, while the thin light [pink] lines show the first (second) order expansion of wave function 
(potential). The dotted line gives the third order expansion of the potential . The dashed light 
[cyan] curve gives analytical approximate forms based on constraints. In the insert, the [red] dashed 
curve gives the potential and the density of the wave-function as offseted by the eigenvalue K. 

(x + 0(x 2 ), with ( = 0'(O), and the first three terms in the expansion of the potential ensues 



We can also identify U = — < x >. 

The numerical solution of Eq. ([3]) is obtained using iteration, similar to the solution 
of the full SP equation. For a given potential U(x), the Schrodinger part is discretized, 
with uniform grid-spacing Ax according to the finite difference method, and the eigenvalues 
and eigenvectors are determined with a banded eigenvalue solver.— K equals the minimal 
eigenvalue and the its normalized eigenvector gives the wave function (p. Next, using this 
result as input, the potential U is updated until self-consistency is reached. Potential mixing 
secures convergence. To improve accuracy and simplify extraction of parameters, I use 
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asymptotic boundary conditions (abc): Since, the wave function falls of exponentially, a hard 
wall boundary condition at x = L for some large cutoff length L is commonly used; however, 
since only a single wave function is retained, the condition <fi(L + Ax) = exp(—\/KAx)(J)(L) 
can be adopted with Ax being the grid spacing. K is updated alongside the potential U in 
the iterative loop. Unlike the hard wall condition, abc guarantees asymptotic behavior at 
the boundary and </>(L) can be used to obtain A^. 



TABLE I: Convergence study of numerical solution. L is the length of the unit cell, Ax is the grid 
spacing. The table give (1 - K(L, Ax)/K com/ ) ■ 1(T 6 where I set K conv = if (50, 1/400). 



L\h 


1/10 


1/50 


1/100 


1/200 


1/300 


1/400 


30 


2973 


442.1 


196.2 


79.25 


41.12 


22.22 


40 


2414 


330.8 


140.6 


51.46 


22.60 


8.329 


50 


2080 


264.1 


107.3 


34.79 


11.48 






Table [I] shows the result of the convergence study. A grid spacing of Ax = 1/400 and a 
length of L = 50 converges K within 10~ 5 . 



TABLE II: Parameters extracted from the numerical solution (described in text), and correspond- 
ing relations to results for specific parameters, c and d are parameters of the constrained-based 
approximate fit. 



Relation 


Parameter 


Value 


<P'(o) = Vac 


c 


0.70708 


Ho) = u'(oh/x 


U'(0) 


1 


z = Xx 


< x > 


2.2543 


E = ~/K 


K 


-0.25902 


g = 2V-K 


9 


1.0179 


V(0) = jU 


U 


- < x > 


Ay = jAu 


Au 


-5.2444 


A,p = Vx 


A^ 


2.3310 


C ~ i 1 + U g) Au-Uo 


c 


0.98957 




d 


1.4256 
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Figured] displays the general, system- independent wave function and potential U, while 
table [Til summarizes key parameters. There is only a single bound state; its eigenvalue is K. 
The value of Uq gives a general relation between between the charge at the interface a and 
band banding at x = 0: V(0) = U (h 2 e/2ma) a 2 / 3 . 

Modelling of semiconductor surfaces and interfaces can benefit from analytical approxima- 
tive expressions of the wave function <p app and potential U app . I here present such expressions, 
which are based on constraints stemming from the numerical solution and physical principles. 
For small x, I choose the conditions U app (x) = Uo + x + 0(x 2 ), and apP (^) = C x + 0( x2 )) 
where ( = <fi'(0), while for large x the approximative expressions should have the exact 
asymptote. Moreover the wave function (j)a, pp (x) should be normalized, which follows from 
charge neutrality. The nodeless shape of the wave function and potential in Fig. [T] motivates 
simple expressions: 

77 ( \ - A U U ° e ~ 9X 

app[X> U + (Au-U )e-^' {b) 
CA xe~ gx / 2 

which obeys the specified conditions if c = (1 + l/U^g) [Au/{Au — Uo)} and d is adjusted 
to normalize the wave function. Table [TTJ lists the determined values of c and d. In Fig. [T] 
the dashed light curves give the approximative solution, which differ from the numerical by 
less than the width of the curves. Fig. [2] details the relative difference: (U — U app ) /U &vv = 
AU/U app and (<p — v 9 a PP )/v :; 'app = A<^/y9 ap p, which is less than 1.6 % for (p and 1 % for 
U. The insert shows absolute differences. The tiny discrepancy between the approximative 
and the numerical solution makes the analytical expressions sufficient for most modelling 
purposes. It also shows that constrained-based strategies can lead to excellent approximative 
expressions. 

The general solution has a limited range of validity. However, it can still serve as an 
approximate account of band bending at surfaces and interfaces capturing essential trends 
and as a building block in simple models. To make a specific test of its robustness, I 
consider an interface charge of o = 10 13 cm -2 , and bulk that is either undoped or doped to 
Pd = 10 19 cm~ 3 , with other parameters as in GaN.— 

Fig. [3] displays the result of the robustness test, which show that for this large charge, the 
general solution is a much better approximation for the potential when bulk is doped, that 
when it is undoped. This result can be understood as follows: As the first excited state gets 




FIG. 2: Relative difference between the constraint-based approximation and numerical solution for 
potential U and wave function 92. The dark [black] curve gives (U — C/ app )/C/ a pp = AC//C/ ap p, while 
the light [cyan] curve gives (ip — v?app)/v?app = A</?/</? ap p. The insert shows absolute differences. 

occupied, it localizes and significantly increases the effective screening length. With doping, 
the excited states instead tend to delocalize, and their contribution to the negative charge 
density partly cancels with the background doping. The complete screening of the interface 
charge by the ground state is therefore equivalent with the approximation that the charge 
density of the excited states cancels with that of the bulk donors, which is impossible for zero 
or tiny donor density. The ground state energy does not typically agree well with that of the 
general solution. The insert shows that for a fairly small interface charge a = 10 12 /cm 2 and 
low temperature, the general solution agrees well with the solution of the full SP equation 
yields virtually identical potentials. This agreement also verifies the accuracy and robustness 
of the numerical solver used and presented in Ref. 

Finally, Fig. @] illustrates the usefulness of the approximative analytical form U app given 
in Eq. ([6]). It shows the potential profile for a system of a two nanometer wide A1N barrier 
sandwiched between GaN cladding layers at different bias Vbias with bulk doping p d as before. 
The full curves give the results of a simple model based on the general solution, while the 
dashed, the result of a full SP calculations using a step-function for the effective Fermi level. 
Such profiles are often displayed only for zero bias, perhaps because of the computational 
complication introduced by a nonconstant Fermi level. To the left of the barrier, a depletion 
layer forms,— ^ which in the simple model is accounted for by a homogeneous charge density 
equalling the donor density pd with a depletion length of L dep . The analytical approximation 
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FIG. 3: Assessment of the general solution beyond its range of validity. The full curves give the 
potential, while the dashed indicate the wave-functions squared (in arbitrary units offseted by their 
energy.) The upper [black] curves give the results obtained from the approximative form of the 
general solution using a sheet density of a = 10 13 cm -2 , and m = 0.2m e . The lower [gray] curves, 
the result obtained using a full Schrodinger-Poisson (SP) solver as described in Ref. Q for the same 
parameters. The (middle [pink]), includes bulk doping of pd = 10 19 cm~ 3 . The insert compares the 
full SP result with that based on the general solution for a = 10 12 cm -2 at zero temperature. 

for the potential profile 7£/ app (a;/A) describes the the inversion layer at the right, with an 
energy 7 and length scale A that depends on the charge of the 2DEG c^deg- The bias Vbias 
over the structure determines this charge: 



where <7 po i is the charge stemming from the spontaneous and piezoelectric effects. Charge 
neutrality gives pd-^dep = o"2Deg- The full potential profile follows from simple electrostatics. 
The result using this model agrees well with that of the SP calculation, which shows that 
the general, system- independent, solution can provide a quick and fairly accurate account 
of how the polarization in AIN/GaN heterostructures influence the potential profiles. 

In summary, the SP equation for a charged hard wall has, in the limit of complete 
screening by the ground-state, been expressed as a general, dimensionless equation. It leads 
to a simple relation between the charge at an interface and the conduction band bending. 
The approximative analytical expressions based on constraints stemming from the exact 
solution provide a convenient tool for obtaining the potential profile and charge density of 
heterostructures with large interface charges, such as for AIN/GaN structures. This could 
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FIG. 4: Potential profile for A1N barrier sandwiched between doped (pd = 10 19 cm~ 3 ) GaN cladding 
layers obtained from a SP calculation (dashed [cyan] curve) with a step-function for the effective 
Fermi level and using an approximative analytical form (full [black] curve) for different biases. 

aid the design of intersubband devices in these materials. 
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